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Recently, there have been significant new insights concerning conditions 
under which closed systems equilibrate locally. The question if subsystems 
thermaUze — if the equilibrium state is independent of the initial state — is 
however much harder to answer in general. Here, we consider a setting in 
which thermalization can be addressed: A quantum quench under a Hamil- 
tonian whose spectrum is fixed and basis is drawn from the Haar measure. If 
the Fourier transform of the spectral density is small, almost all bases lead to 
local equilibration to the thermal state with infinite temperature. This allows 
us to show that, under almost all Hamiltonians that are unitarily equivalent 
to a local Hamiltonian, it takes an algebraically small time for subsystems 
to thermalize. 
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1 Introduction 

Consider a large closed system suddenly taken out of equilibrium and evolving 
unitarily in time. Now suppose access to the system is limited and only a small 
subsystem may be observed. If this subsystem equilibrates and information about 
the initial state is lost, principles of statistical mechanics emerge locally: Given 
sufficient time to equilibrate, one would observe a generic state — e.g., a (gen- 
eralized) Gibbs state — which does not depend on the fine-grained information 
contained in the initial state anymore but, possibly, instead on macroscopic ob- 
servables such as, e.g., the mean energy. The system under consideration is closed 
and the dynamics are unitary. Ignorance — the full extent of which one ought to ac- 
knowledge frankly [|T1|, arriving at Jaynes' principle — is replaced by limited access 
(spatially - only subsystems may be observed) to the system. There is no need 
to introduce ignorance - no need to put probabilities by hand. Instead, generic 
ensembles emerge locally from unitary evolution. This link, along with the re- 
cent experimental availability of such quantum quench settings All [3l IH |51 [6l |3, 
has spurred a bit of a renaissance of ideas dating all the way back to von Neu- 
mann's quantum ergodic theorem |[8l. An incomplete list of theoretical studies of 
these questions include the eigenstate thermalization hypothesis [l9l[T0l, solvable 
systems [[III [El [HI [111 [151 [l6l, conformal field theoretical settings and nu- 
merical work [[181 . Conditions guaranteeing and ruling out thermalization have 
also been given recently [[T9l[20l[2T]| . 

For several reasons, the question of equilibration after a quench to a random 
Hamiltonian has attracted quite some attention lately [[22l [23l [24ll : A random 
Hamiltonian may serve as a model for a sufficiently complex system (e.g., certain 
properties of quantum chaotic systems are well described by random Hamiltoni- 
ans [[241 [251 [261 ). In fact, in Ref. [[231 , a relation between complexity (defined 
as the number of one- and two-qubit gates needed to approximate the unitary 



2 



diagonalizing the Hamiltonian) and equilibration properties was suggested. Ad- 
ditionally, results on random Hamiltonians have implications on the time scale on 
which random quantum states may be generated [24J. Furthermore, to the best of 
our knowledge, no explicit analytic results concerning the equilibrium state and 
equilibration time are known for non-integrable systems. As we will see, the set- 
ting of random Hamiltonians puts us not only in the position to derive explicit 
bounds, but also allows us to make the equilibrium state and the time scale of 
equilibration very explicit. It is also a setting in which usual assumptions on the 
degeneracy of the energy spectrum are not necessary. 

In more technical terms, we are concerned with the following setting. The sys- 
tem is initially in the state \ipo) (below we will also allow for mixed initial states) 
and evolves unitarily according to a Hamiltonian H, g(t) = e~^^^ \tpo) {tpo\e^^^ . 
We distinguish a subsystem S and consider the state on S, fc(t) = tTQ[g{t)], 
where B denotes the rest of the system, which takes the role of a "bath" or 
"environment". Now we ask the question wether there is some canonical time- 
independent state oj which describes the system locally, i.e., wether Qs{t) is close 
to Us := trQ[uj], and if so, for which times. That is, we are concerned with the 
trace distance ||^5(t) — aj^Htr, which quantifies the distinguishability of the two 
states: For any observable O on 5 one has \{0)gg(t) — {0)^^\ < ||^5(t) — 
li'^lltr, where || ■ || denotes the operator norm. A seminal result is the following ll27l 
(see also Ref. [|28]| ). If the energy gaps of the Hamiltonian H are non-degenerate, 
one has 



where is the dimension of the Hilbert space associated with S and tr[ii;^] is the 
purity of the equilibrium state 



where \vn) and are eigenstates and corresponding energies of the Hamiltonian, 
respectively, and tpn = {vn\ipo). Hence, this result identifies the equilibrium state 
and establishes that if its purity is small, the system is locally well described by uj 
for almost all times. Two questions remain. What is the time-scale of equilibra- 
tion, i.e., what can be said for finite times T? Further, when is ujs independent of 
the initial state, i.e., when does the system not only equilibrate but thermalize in 
this sense? Progress towards the equilibration time scale was made in the recent 
Ref. |[29l , in which also the condition on energy gaps was relaxed, arriving at a 
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bound that also involves the purity of the equilibrium state. Picking the Hamilto- 
nian H = Udiag[{En}]W randomly by fixing the spectrum {En} and drawing 
U from the Haar measure, conditions on the energy gaps may also be relaxed and 
equilibration time scales obtained [|22ll23l . The question of initial state indepen- 
dence, i.e., thermalization remains open, however. For the randomized setting just 
described, we will show that thermalization can be addressed and that the system 
thermalizes to the maximally mixed state = l/dg - the Gibbs state with in- 
finite temperature. If the spectrum {En} is that of a local Hamiltonian, we will 
show on which time scales to expect thermalization. 



2 Thermalization under random Hamiltonians 

We start with some notation. The Hilbert space under consideration is 

"H = span{|n) = \ni ■ ■ ■ un) ni = 1, . . . , di} = S ®B 

and we denote ds = dim[iS], d^ = dim[i3], d = dim['H] = dgd^. This is the 
Hilbert space of a collection of N subsystems, each with dimension di - e.g., 
spin- 1/2 particles (di = 2). Until we consider local Hamiltonians below, all 
considerations are independent of the geometry of the modelled system. That is, 
we need not assume a specific arrangement of the N subsystems or sites (such as, 
e.g., spins on some sort of lattice). Let H a Hamiltonian, 

^ = ?7($^^n|n)(n|)?7t, (3) 

n 

g{t) a time-evolved state, 

m = e-'*^^oe'*^, (4) 

and denote Qs{t) = trB[^(t)]. We do not restrict ourselves to pure initial states 
but allow ^0 to be mixed. We will pick Hamiltonians as in Eq. ^ at random by 
fixing the spectrum {En} and sampling U from the Haar measure. We denote 

E[.] = j»dfi{U), (5) 

where yu is the Haar measure. We set out to bound the expected trace distance of 
the local time-evolved state Qs{t) to the maximally mixed state. If this distance 
is small, the system is locally well described by the thermal state with infinite 



4 



temperature. In Ref. [l24ll it was shown that for initial pure product states, = 
iV's) ® \^b) the expected purity of Qs{t), E [tr5[^|(t)]] , is given by JSOj 

j_ A{d-5+i)<p Ami , , n^\twm] _ m)? \ 

l + d (d + 3)(rf2-l)l A Ad" 2d d^ J' 

(6) 

where 6 = d^ + ds and 

k 

is the Fourier transform of the spectral density. For general states g, the quantity 
(pg{t) = tr[e'*^^] is also known as characteristic function (it is a positive definite 
function and hence, due to Bochner's theorem, the characteristic function of a 
random variable), taking centre stage in many proofs of quantum central limit 
theorems (see, e.g., Refs. ll36l l37l ). and is the Loschmidt echo of \ip). 

Statistics and equilibration time of the latter was studied recently in [|3T1l . 

As \4>{t)\ < 1, it follows for separable initial states = Tl,nPnQ~s'® Qb^ that 

and we may formulate the following direct consequence of which should be 
compared to Refs. [|22ll23l , in which a similar bound was given for the expected 
distance of Qs{t) to the (there unknown) equilibrium state in Eq. Q. 

Corollary 1 Let {En\ be given, a separable state, and H as in Eq. Then 

^^\Qs{t) - Is/dsWtr] < ^s^WW^Wb- (8) 

Hence, if the right hand side is small, we expect the system to be close to the 
maximally mixed state - the thermal state with infinite temperature. 

We now turn to investigating (pit) for two classes of systems for which the 
above bound can be made explicit and also the thermalization time scale be given 
explicitly: Solvable systems and general local Hamiltonians. The latter constitutes 
the main result of this work while the former illustrates — at the hand of a rather 
simple proof — what to expect for more general systems. Due to it being rather 
technical and long, the proof for general local Hamiltonians may be found in the 
Appendix. 
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3 Time scale of thermalization 



We write (■) = tr[^ ■] and 

a^ = {{H-{H)n -^ = ^, (9) 

where we note that for a large class of Hamiltonians a is simply a constant ll33l . 
Further, we denote the time average of the trace distance as 

A(T) = ^ ^ Usit) - Is/dsWtr. (10) 

We will now give explicit bounds on E[A(T)] in terms of a and the system size, 
enabling us to extract the time scale of thermalization. Suppose we obtain the 
bound E[A(T)] < c. We may then also give a bound on the probability that 
A(T) < yc: From Markov's inequality, we have for all ?/ > that 

F[A{T)<yc] >l-l/y. (11) 

Further, the fraction of times in [0, T] for which — l^/fi^Htr < l/a^ is at 

least 1 — a;A(r). Combining all the above, a bound E[A(T)] < c hence enables 
us to conclude that with probability at least 1 — 1/y, the fraction of times in [0, T] 
for which Wgsit) — ls/ds\\tv < 1/xis at least 1 —xyc. Hence, choosing x,y 1 
and supposing c is such that xyc ^ 1, we have that for almost all U and almost 
all t G [0, T], the subsystem is close to the maximally mixed state. We proceed 
by giving bounds for solvable system and then show that similar bounds hold for 
spectra of local Hamiltonians. 



3.1 Solvable systems 

Assume that the spectrum takes the form 



N 



= ^ efcTifc, Uk = 0, 1. (12) 

fc=l 

This includes, e.g., spin chains solvable via a Jordan-Wigner transformation, i.e., 
in particular all Hamiltonians of the form 

N~l N 
i=l a,fS=x,y i=l 
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i.e., all XF-type models in transverse fields. Straightforward calculations show 
that 



4 

k 



and 



|0(t)| = nicos(efct/2)|. (15) 



k 



Hence, denoting emax = max^ |efc|, we have for |t|emax < 27r that 

k 

which implies that there is an absolute constant such that for Temax < Svr, we 
have 

E[A(r)]^<M.(^ + i). (17) 

In particular, 

E[A(rr<M.(i^.i-) (18) 

for Temax = X*""^/^ and all < e < 1/2. Hence, for sufficiently large system 
size and almost all Hamiltonians as above, the subsystem will spend most of its 
time in [0,T], T ~ N''"^^'^, close to the maximally mixed state - under almost 
all Hamiltonians with a spectrum as in Eq. (fT2l) . subsystems thermalize in a time 
T ~ N"'^/'^. We will now see that a similar result holds for spectra of local 
Hamiltonians. 



3.2 Local Hamiltonians 

We now consider local Hamiltonians on a /^-dimensional cubic lattice, which we 
denote as the collection of sites £ = {1, . . . , , i.e., N = M^. We equip the 

lattice with a distance dist{i,j), which we take as the shortest path connecting i 
and j. For open boundary conditions, this is simply dist{i,j) = J2f=i I'^s — jsl- 
The Hamiltonian is assumed to be local in the sense that 

H = J2h^^ (19) 

iec 

where hi acts non-trivially only on {j E C \ dist{i,j) < R} for some constant 
< -R G N. Further, we assume that \\hi\\ < h for all i and some constant h and 
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that the lattice is sufficiently large such that < < M^/^. Hamiltonians as 

just described, we will simply call local Hamiltonians. A proof of the following 
theorem can be found in the appendix. The proof is based on ideas (which we 
generalized to account for finite system sizes N) employed in proofs of quantum 
central limit theorems Il36ll37]| . 

Theorem 1 Let H a local Hamiltonian. Then there exist constants and bo, 
only depending on the interaction radius R, the interaction strength h and the 
dimension of the lattice D, such that for 

T = aoa2iVi/(5^)-i/2 (20) 



we have 



ljym\<^^- (21) 



This together with ([8]) allows us to formulate our main result. 

Corollary 2 Let {En} be the spectrum of a local Hamiltonian. Let H as in Eq. 

and Qq a separable state. Then there are constants oq (^'^d bo, only depending 
on the interaction radius R, the interaction strength h and the dimension of the 
lattice D, such that for 

T = aoa'iV'/^'^)-'/' (22) 



we have 



( 



yNi/{5D) + 



Hence, under almost all Hamiltonians that are unitarily equivalent to a local Hamil- 
tonian, the subsystem will, for sufficiently large system size, spend most of its time 
in [0, T], T ~ N^/i^D)"'^/'^ ^ close to the maximally mixed state - under almost all 
Hamiltonians that are unitarily equivalent to a local Hamiltonian, subsystems ther- 
malize in a time T ~ N'^/i^^)-'^/'^ _ More precisely: Let e > and T as in Eq. 
(|22l) . Then, with probability at least 

1 _ n^^-jd)/\ (23) 

the fraction of times in [0, T] for which 

\\gsit)-ls/ds\\tr<N^'-^^^^ (24) 
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is at least 



1 




1 + a-3 iVi/(5i3) 



) 



(25) 



Maybe surprisingly, the size of the subsystem does not need to be constant. The 
above bounds allow for subsystems whose size increases logarithmically with the 
system size A^. 

4 Discussion 

We have shown on which time scales subsystems thermalize under unitary dy- 
namics generated by randomized local Hamiltonians. In this setting, usual as- 
sumptions on the degeneracy of energy gaps are not necessary and the equilibrium 
state — here the maximally mixed state — and bounds can be given explicitly. The 
only remaining assumption is that the energy variance is lower bounded by the 
system size, which is also necessary for asymptotic normality Il36l 1371 and ful- 
filled for a large class of Hamiltonians ll33l . As the system locally equilibrates to 
the maximally mixed state — the thermal state with infinite temperature — the ques- 
tion arises under which conditions a thermal state with finite temperature might 
emerge. It would be interesting to see wether restricting the unitaries to the ones 
preserving the mean energy of the initial state lead to a finite temperature Gibbs 
state. 
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A Proof of Theorem 1 
A.l Preliminaries 

We denote (■) = tr[-]/(i. For subsets ^ C £, we write 



and omit the index for A = C. We will first prove the following lemma and then 
make the partition explicit in section [A4l which proves theorem 1. 

Lemma 1 Let C = A U C, A = U^^^An a partition of the lattice such that 
dist{An, Am) > 2Rfor n ^ m. Then there are constants Oq and bo, only depend- 
ing on the interaction radius R, the interaction strength h and the dimension of 
the lattice D, such that for 



T < 



A\ maxn |.4„|i/2 



(26) 
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we have 



i / dt\<p{t)\<bo{\C\T' + ^). (27) 

If the lemma holds for (hi) = then, by noting that 

= |0(t)e-'*S»"<^'>|, (28) 
it holds also for (hi) ^ 0, i.e., we may assume that (hi) = 0. Now, 



(29) 



We proceed by bounding 6i{t) and 62{t). The proof techniques are similar to the 
ones in Ref. [[361 (see also Ref. [[371 ), we however improve on some of the bounds. 
Several times we will make use of 

max V 1 < 1 + c^r-^ =: /3r, (30) 

dist{i,j)<r 

the constant cd depending only on the dimension of the lattice. 
A.2 Bound on 



Denote He = YlnLi ^n, i-C-^ we slightly abuse notation by using i e C and 
1, . . . , |C| interchangeably. Writing Hm = H — X]™=i ^n, Hq = H, we have 

\c\ 



n 



m=l 

m=l 
\C\ 



m=l 

where 



XUt) = XmiO) + t^{0) + fdrTds ^{s), (32) 

Jo Jo 
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-^m(O) = 0, ^#-(0) = ihm, and the last term is equal to 



'o Jo 



2 \„-isH„ 



I.e., 



where (recall that /ij and hj commute for dist{i,j) > 2R) 

\\[hra,Hrn]\\< WK. < (32R. 



dist{i,m)<2R 



Now, Hrn = TjiaXm ^« some C C, i.e., writing 

dist{i,m)<2R dist{i,m)>2R 



we have (/ime = 0, i.e. 



fds\\Yj\<h^P,n\t\. 



Therefore, < |C|^ (1 + S/^si?) 



A.3 Bound on 62 

Writing 



we have 



n=l 
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Taylor's theorem implies 



.W-l|<(^5f =:^, (36) 



g(<x2-a3^)t2/2gEn=l('^^tV2+log[<^n(i)]) _ ^ 



i.e., if cr^t^ < 2, we may take a logarithm and, noting that cr'^ = J2n=i write 
Where F(t) = Ell /n(t), 

/nW = l4^ + log[0n(t)]| 

< \<Pn{t) - 1 + + |</.„(t) - 1 - log[0„(t)]| 

where we used Taylor's theorem, Cauchy-Schwarz, and that for all 1^; — 1| < 1 
one has 2\z - log(l + z)\ < \z\^''^. Now, {H^) < and al < C2IAI (see 

Section |A3] below) . Hence, 

- 4^ + /nW < v^lXI^^/^f - (37) 
For all T > 0, we have 

|C|T^ + ^>(S^)'''> (38) 
i.e., for 8IC2ICI > 4cr^, the lemma is trivially true: 

i / dt|0(t)| < 1 < i^Y^' < i3c,Y/'\C\T^ + ^ 
Jo 

for all T > 0. We may hence assume that 81c2|C| < ia^. As - cr^| < Sc2\C\ 
(see Section lA3] below), we hence have 23(T^ < 27a j^. Now, for x > 

1^1 - max! 1 (39) 

implies 

< C2|A|t^ < ri^^^S^ < clx^ (40) 
15 



and 

^'^^-1'--^^^^ + Fit) < + V^,xa% - 4^ 

/23 

6 6 



Hence, setting x = min{yT7c|, 22/(6^0^)} implies cr^t^ < 2 and 

^l-^-^J-^^^^ + < (41) 
Thus, for t as in Eq. (l39l) , we have 

A.4 Putting it together 

So far, we have that there are constants oq and bo, only depending on R, h, and D, 
such that for 

we have 

/■oo 

<\C\T' + ^ dt(t2 + l)e- 

^0 



r/36 



which proves the lemma. To prove theorem 1, we choose the sets \An\ as follows. 
Let K = , a = M/K - 2R, and for A; = 0, . . . , A' - 1, define X = 

{1, . . . , M}^(^-i), Ofc = (a + 2R)k, 



Ak = X X {n e ¥i \ ak < n < ak + a} 

= X X {\ak-l],...,[ak + a\}, 
Ck = Xx{nE¥l\ak + a<n< flfc+i} 



(43) 



= X X {\ak + a-l],..., [ofc+ij }, 
C = |jf=o^ Ck- Then, for n^m, 

dist{An, Am) > K+iJ - [ofc + a - 1] + 2 > 2i? + 1. (44) 
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Assuming M^/^ -2R> and M - 2M^/^ > M/2 implies 



and 



^ M3/5+2fl <K<M ' 



2 < M'^-'a = M'^-^M/K -2R) 

< \Ak\ = M^-i(Lafc + aj - K -11+1) 

< M^-\M/K -2R + 2)<^ 



(45) 



(46) 



Further, 



Hence, putting 



\Ck\ = M^-i([afc+iJ - K + a - 1] + 1) 
< 2{R+l)M^-\ 



T 



Oo 



mV 5-o/2 _2 ^ 



|^|max„|yt„|i/2' 



we have, redefining bo. 



4 /dt|0(t)|<M-V5 + ^^^^. 



(47) 



(48) 



Finally, the conditions M - 2M^/^ > M/2 and M^/^ - 2i? > are implied 

by < 43/27? < which also implies K > 1, | A| > 1, and \Ck\ > 1. 



A.5 Some bounds 

In this section we will derive the bounds cr^ < C2|^n|, {H^) < C4|^„p, and jcr^ — 
cr^il < 3c2|C|. We recall that {hi) = OforalH,that (hihj) = for dist{i, j) > 2R, 
and the bound in Eq. (|30] ). We have 

dist(i,j)<2R dist{i,j)<2R 

Similarly, 

a^-a\= {Hi) - 2{Hj,Hc) < C2\C\ - 2{H^Hc), (50) 
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where 



ieA,jeC 

dist{i,j)<2R 



(51) 



i.e., |cr^ — cr^l < 3c2|C|. To show that < C4|^„p, we define for given i G An 

Hn=Y. ^^^= Yl + Yl =• + (52) 

j&An jeAn jeAn 

dist(i,j)<2R dist{i,j)>2R 

for which we have [hi, Yi\ = and (hiYj^) = 0, i.e., 
ieA„ 



ill/) 



(53) 



where we used the cyclic property of the trace. Now define 

5Z ^.■= Yl ^^■+ Y k=■R^ + Z^ (54) 

j&An j&An jeAn 

dist{i,j)>2R 2R<dist(i,j)<AR dist{i,j)>4R 

for which we have [hi, Zi] = 0, [Xi, Zi] = 0, and (hiXiZf) = 0, i.e., due to the 
cyclic property of the trace. 



J2 {kXiY^) = Y ChiMRl + '2RiZi)) 

i£An ieAn 

</.5^ (||X,||||i?,f + 2||X,||||^,||||Z,||). 

i&An 



(55) 



Noting that 

11^.11 < hf32R, IIF.II < h[An[, [[Z4 < h[An[, m\[ < h(3^R, (56) 



we finally have 

{H^} < h' [{3P2RPln + (3lR)\An[ + mRf^m + Sf3l^)\An['] 
< h" [3P2rPIr + PIr + 6l32Rl3m + 3(3^^^] [An[' =: C4IA 



(57) 
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